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1 Abstract 

2 

3 We use observations from two aircraft during the ICARTT campaign over the eastern United States 

4 and North Atlantic during summer 2004, interpreted with a global 3-D model of tropospheric chemistry 

5 (GEOS-Chem) to test current understanding of the regional sources, chemical evolution, and export of 

6 nitrogen oxides. The boundary layer NO x data provide top-down verification of a 50% decrease in power 

7 plant and industry NO x emissions over the eastern United States between 1999 and 2004. Observed 8-12 

8 km NO x concentrations in ICARTT were 0.55 ± 0.36 ppbv, much larger than in previous United States 

9 aircraft campaigns (ELCHEM, SUCCESS, SONEX). We show that regional lightning was the dominant 

10 source of this NO x and increased upper tropospheric ozone by 10 ppbv. Simulating the ICARTT upper 

1 1 tropospheric NO x observations with GEOS-Chem require a factor of 4 increase in the model NO x yield 

12 per flash (to 500 mol/flash). Observed OH concentrations were a factor of 2 lower than can be explained 

13 from current photochemical models, and if correct would imply a broader lightning influence in the upper 

14 troposphere than presently thought.. An NO y -CO correlation analysis of the fraction / of North American 

15 NO x emissions vented to the free troposphere as NO y (sum of NO x and its oxidation products PAN and 

16 HN0 3 ) s shows observed /=16±10% and modeled /=14±8%, consistent with previous studies. Export to 

17 the lower free troposphere is mostly UNO;, but at higher altitudes is mostly PAN. The model successfully 

18 simulates NO y export efficiency and speciation, supporting previous model estimates of a large U.S. 

19 contribution to tropospheric ozone through NO x and PAN export. 
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1 1. Introduction 

2 Quantifying the sources and fate of nitrogen oxides (NO x = NO + N0 2 ) over northern mid- 

3 latitudes continents is critical for assessing anthropogenic influence on global tropospheric ozone 

4 [j Pickering et al., 1992; Jacob et al., 1993; Thompson et al., 1994; Li et al., 2004]. The 

5 International Consortium on Atmospheric Transport and Transformation (ICARTT) aircraft study 

6 [Singh et al., 2006a; Fehsenfeld et al., 2006], which took place in July - August 2004 over the 

7 eastern United States and the North Atlantic, provides an opportunity for this purpose. We present 

8 here a global 3-D model analysis of ICARTT observations for NO x , its chemical reservoirs, and 

9 related species including hydrogen oxide (HO x ) radicals to quantify continental NO x sources from 

10 combustion and lightning, determine the chemical fate of NO x in the United States boundary 

1 1 layer and in North American outflow, and examine the implications for ozone. 

12 Ozone production in the troposphere is principally limited by the supply of NO x [Chameides, 

13 1992], Fossil fuel combustion accounts for over half of the global NO x source [Intergovernmental 

14 Panel on Climate Change (IPCC), 2001]. The United States has been actively reducing its 

15 summertime NO x emissions since 1998 to decrease ozone smog. The Environmental Protection 

16 Agency (EPA) NO x State Implementation Plan (SIP) Call mandated that 22 eastern states meet 

17 state-specific total NO x emissions reductions by 2003 (Phase 1) and further reduction by 2007 

18 (Phase 2). By 2003, all 22 states had reduced NO x point source emissions to their Phase 1 levels. 

19 Frost et al., [2006] determined from stack sampling that power plant emissions of NO x decreased 

20 50% between 1999 and 2003. NO x levels in the U.S. will likely continue to drop. In March 2005, 

2 1 the EPA issued the Clean Air Interstate Rule, which will, when fully implemented, permanently 

22 reduce NO x emissions to 60% of 2003 levels in 25 eastern states. The ICARTT observations 

23 offer an opportunity to check on these emission reductions. 

24 Oxidation of NO x to HN0 3 , peroxyacetylnitrate (PAN), and other minor products takes place 

25 on the order of hours, so that rapid ozone production is confined to the continental boundary layer 

26 (CBL), seemingly limiting its affect on global ozone. However, the dependence of ozone 
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1 production on NO x is highly nonlinear; the ozone production efficiency per unit NO x consumed 

2 (OPE) increases rapidly as the NO x concentration decreases [Liu et al., 1987]. This means that a 

3 small fraction of emitted NO x exported to the free troposphere by frontal lifting, deep convection, 

4 or boundary layer venting could lead to significant ozone production in the free troposphere over 

5 the continent or downwind [Jacob et al., 1993; Thompson et al., 1994], Similarly, PAN (which is 

6 thermally unstable and not water-soluble) can be vented from the boundary layer and transported 

7 on a global scale at cold temperatures, eventually decomposing to release NO x as air masses 

8 subside. Quantifying the sources, chemical evolution, and export of anthropogenic NO x (and 

9 PAN) is thus critical to understanding the North American contribution to the global ozone 

10 budget. This contribution is important from the perspective of ozone as a greenhouse gas 

1 1 [Mickley et al., 2004] and for intercontinental transport of ozone pollution [Jacob et al., 1999; 

12 Holloway et al., 2003]. 

13 Early Eulerian models found that the fraction /of NO x emitted in the United States that is 

14 exported out of the CBL as NO y (sum of NO x and its oxidation products) ranges from 25% in 

15 summer to 35% in winter [Horowitz et al., 1998; Liang et al., 1998]. Subsequent Lagrangian 

16 analyses using NO y -CO correlations measured from aircraft in free tropospheric outflow (2-6 

17 km), over the North Atlantic in September (NARE ’97) seemed to contradict these model results, 

18 with estimates of / ranging from only 3% [Stoh l et al., 2002] to 9 ± 14 % [Parrish et al., 2004], 

19 Li et al. [2004] reconciled the Eulerian and Lagrangian approaches by pointing out that the early 

20 Eulerian models had insufficient UNO; scavenging, while the Lagrangian models underestimated 

21 background CO. They derived a consistent value /= 17-20% by both approaches for the 

22 NARE’97 period. The ICARTT study offers far more geographical coverage and chemical 

23 information in the boundary layer and the free troposphere than previous studies, enabling better 

24 constraints on the estimates of anthropogenic export and associated NO y speciation, as well as the 

25 underlying source and chemical factors. 
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1 In addition to convectively lofted NO x , a highly uncertain source of NO x to the upper 

2 troposphere is from lightning. Global lightning source estimates range from 1 to 12 Tg yr' 1 [ Price 

3 et al., 1997; Nesbitt et al., 2000], with the most recent estimates in the range 1-6 Tg N yr 1 

4 [ Boersma et al., 2005; Beirle et al., 2006; Martin et al., 2006b], Past studies disagree on the 

5 relative importance of lightning versus convective injection in supplying upper tropospheric NO x 

6 [Jaegle et al., 1998; Levy et al., 1999; Grewe et al., 1999; Li et al, 2005]. Here we use the 

7 ICARTT data to examine the contribution of each. 

8 2. The ICARTT Study 

9 ICARTT took place over eastern North America and the North Atlantic in July-August 2004. 

10 A major objective was to quantify North American sources and outflow of pollutants and 

1 1 climatically important species. Two principal components directed at that objective were the 

12 NOAA New England Air Quality Study/Intercontinental Transport and Chemical Transformation 

13 (NEAQS-ITCT 2004) and the NASA Intercontinental Transport Experiment - North America, 

14 Phase A (INTEX-A). 

15 The NOAA NEAQS- ITCT 2004 campaign [ Fehsenfeld et al., 2006] took place July 3 - 

16 August 15, 2004 over the NW Atlantic and the NE United States out of Portsmouth, New 

17 Hampshire (Figure 1). It used a WP-3D aircraft (ceiling ~6 km). The NASA INTEX-A campaign 

18 [ Singh et al., 2006a] took place June 29 - August 14 over the central and eastern United States, 

19 and the North Atlantic, from bases at Edwards (California), St. Louis (Missouri), and Portsmouth 

20 (New Hampshire). It used a DC-8 aircraft (ceiling ~12 km) with extensive vertical profiling and 

2 1 boundary layer mapping at 1 kft. (Figure 1) The WP-3D remained close to New England 

22 throughout the mission, focusing on emissions verification and chemical transformation of major 

23 urban pollution plumes. Detailed descriptions of aircraft payloads and measurement techniques 

24 are in Fehsenfeld et al., [2006] and Singh et al., [2006a] for the WP-3D and DC-8 respectively. 

25 We make use here of the following 1 -minute average measurements (accuracies given in 

26 parentheses for DC-8, WP-3D respectively): ozone (5%, 3%); CO (5%, 2. 5%); CH 2 0 (10%, DC- 
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1 8 NCAR measurement only),; HN 04 ( 15 %, DC -8 only); H 2 O 2 ( 20 %, DC -8 only); NO (15%, 5%) 

2 DC -8 Pennsylvania State U. measurement only); N0 2 (10%, 8 %); HN0 3 (15%, 15%); 

3 PAN(15%,15%), OH(15%, DC -8 only), and H0 2 (15%, DC -8 only). We estimate the 

4 concentration of NO y in free tropospheric export from the DC -8 data as the sum of NO x , PAN, 

5 and HN0 3 . Nitrate aerosol, organic nitrates other than PAN, and HNO 4 were also measured 

6 aboard the DC -8 [Cohen et al. , 2000; Dibb et al., 2003; Huey et al., 2004], but accounted for less 

7 than 10% of NO y in the free troposphere and the data sets were sparse. 

8 Ventilation of the eastern United States in summer is primarily driven by cyclones 

9 tracking eastward typically every 5 days in the 45° - 55°N band [Li et al., 2005]. During ICARTT 

10 a persistent trough along the east coast led to cyclones extending further south [Fuelberg et al., 

1 1 2006], and lack of stagnant surface high pressure led to a record low number of air quality 

12 violations [Thompson et al., 2006]. Parrish et al. [2004] found that outflow from fair weather 

13 cumulus rising to 3-4 km provides an additional ventilation mechanism. Yet another mechanism, 

14 particularly in the south, is deep convection from surface heating [Li et al., 2005]; this mechanism 

15 is often associated with lightning activity. 

16 3. Model Description 

17 We simulate the ICARTT observations with the GEOS-Chem global 3-D model of 

18 tropospheric chemistry (version 7.02; http://www.as.harvard.edu/chemistry/trop/geos/) driven by 

19 assimilated meteorological observations from the Goddard Earth Observing System (GEOS-4) of 

20 the NASA Global Modeling and Assimilation Office (GMAO). The model is applied to a global 

21 simulation of ozone-NO x -hydrocarbon-aerosol chemistry with 120 species simulated explicitly. A 

22 general description of GEOS-Chem is given by Bey et al. [2001] and a description of the coupled 

23 oxidant-aerosol simulation as used here is given by Park et al. [2004], Partitioning of total nitric 

24 acid between the gas and aerosol phases is calculated using the MARS-A thermodynamic 

25 equilibrium model [Binkowski and Roselle, 2003]. Emissions in the model are as described by 

26 Park et al. [2004] unless specified otherwise. 
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1 Meteorological fields in the GEOS-4 data have a temporal resolution of 6 hours (3 hours 

2 for surface variables and mixing depths) and a horizontal resolution of 1° latitude by 1.25° 

3 longitude, with 55 vertical sigma levels between the surface and 0.1 hPa (including about 16 in 

4 the troposphere and 5 in the boundary layer up to 2 km). For input to GEOS-Chem we degrade 

5 here the horizontal resolution to 2° latitude by 2.5° longitude. Mean afternoon boundary layer 

6 heights over the eastern United States in the GEOS-4 data for the ICARTT period are 1 100 ± 400 

7 m AGL. Tropopause heights are 13.3 ± 2.2 km. The cross-tropopause ozone flux is specified 

8 globally with the Synoz method [ McLinden et al., 2000] while the NO y flux is calculated from 

9 N 2 0 oxidation in the model stratosphere [Bey et al., 2001]. Global net cross-tropopause fluxes of 

10 ozone and NO y are 495 Tg O 3 yr ' 1 and 2 Tg N yr' 1 , respectively. The model wet deposition 

1 1 scheme [Liu et al., 2001] includes contributions from scavenging in convective updrafts, and 

12 rainout and washout from convective anvils and large-scale precipitation, and it allows for return 

13 to the atmosphere following evaporation. Soluble gases are taken up by liquid water on the basis 

14 of their effective Elenry’s law, by ice on the basis of co-condensation or surface coverage, with 

15 species-specific retention efficiencies when droplets freeze, as described by Mari et al., [2000]. 

16 The simulations are conducted for July-August 2004 and are initialized with a 18-month spin-up 

17 simulation. For comparison with observations, the model is sampled along the aircraft flight 

18 tracks and times. The comparisons exclude fresh pollution plumes as diagnosed by NO x /NO y > 

19 0.4 mol mol " 1 or (if NO y is not available) N0 2 > 4 ppbv and altitude < 3 km; biomass burning 

20 plumes as diagnosed by FICN > 500 pptv or CH 3 CN > 225 pptv; and stratospheric air as 

2 1 diagnosed by ozone/CO >1 .25 mol mol" 1 . These filters exclude 8 %, 22% (fresh pollution plumes); 

22 5%, 8 % (biomass burning plumes), and 6 %, 0% (stratospheric air) of the overall data set for the 

23 DC -8 and WP-3D, respectively. 

24 A major focus of our work is to use the ICARTT observations to test and improve U.S. 

25 NO x emission estimates. We will show simulations with “standard” emissions based on a priori 

26 information from the standard version of the GEOS-Chem model, and “improved” emissions that 
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1 reflect the ICARTT constraints and improved understanding. These emissions for the ICARTT 

2 period (July 1 - August 15, 2004) are summarized in Table 1. Original fossil and biofuel 

3 emissions in the United States are from the EPA 1999 National Emission Inventory (NEI99). 

4 They amount to 0.79 Tg N for July 1 -August 15 and 6.2 Tg N annually, with distribution shown 

5 in Figure 2 (left). Transportation accounts for 35%, industry 17%, power generation 26%, and 

6 other sources 22% (mostly non-road vehicles). As we will see in section 5, NO x concentrations, 

7 are overestimated using the 1999 inventory, in the U.S. boundary layer, consistent with reduction 

8 of NO x emissions from power plants by 50% from 1999 to 2004 driven by the NO x SIP Call 

9 [Frost et al., 2006], These reductions are being added to subsequent inventories. Reduction of the 

10 power plant and industry sources results in the improved inventory of Table 1 and Figure 2 

1 1 (right). The ICARTT data also show that CO emissions in the NEI99 inventory are 30% too high, 

12 as discussed below, and we make this adjustment in the improved inventory. 

13 The global lightning source of NO x in GEOS-Chem is computed with the scheme of 

14 Price and Rind [1992] that relates number of flashes to convective cloud top heights as described 

15 in Wang et al., [1998], We scale the NO x production to 4.7 Tg N yr' 1 globally, resulting from 2.7 

16 x 10 9 flashes, and distribute this source vertically following Pickering et al. [1998] (55-75% 

17 above 8 km, up to 23% in lowest km). The resulting U.S. emissions for the ICARTT period 

18 (original inventory) are 0.067 Tg N for the contiguous United States and coastal waters (130- 

19 70°W, 25-50°N), concentrated over the Southwest and along the Gulf of Mexico. As we will see 

20 in section 6, successful simulation of the ICARTT upper tropospheric NO x observations in 

21 GEOS-Chem requires a factor of 4 increase in this source along with an upward shift to the July 

22 mean tropopause height, and this is included in the improved inventory. 

23 Extensive and persistent boreal forest fires took place in Alaska and NW Canada in 

24 summer 2004; in contrast, there were no significant fires in Siberia. We use the daily biomass 

25 burning inventory of Turquety et al. [2006] for North American fires during ICARTT. This 

26 inventory was constructed by combining daily area burned reports from government agencies and 
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1 hot spots detected from space by the MODIS instrument with estimates of fuel loadings and 

2 emission factors depending on the type of ecosystem burned. This inventory for CO was 

3 evaluated against MOPITT columns as described by Turquety et al. [2006]. North American CO 

4 and NO y fire emissions for July 1 -August 15 are 20 Tg CO, 0.3 Tg N. Short-lived hydrocarbons 

5 emitted from fires drive fast conversion of NO x to PAN, slowing down ozone formation [Jacob et 

6 al., 1992] but this is inadequately represented in the model where the fastest-reacting pyrogenic 

7 hydrocarbon is propene [Hudman et al., 2006]. In the improved emission inventory we release 

8 80% of the biomass burning NO x as PAN. 

9 

10 4. Hydrogen oxide radicals and reservoir species 

1 1 Hydrogen oxide radicals (HO x = OH + peoxy radicals) and their reservoirs (H 2 O 2 , CH 2 0, 

12 HNO 4 ) largely define the photochemical environment for NO x oxidation and ozone production. 

13 Figure 3 shows simulated and observed mean vertical distributions of OH, HO 2 , H 2 O 2 , CH 2 0, and 

14 HNO4 for the ensemble of INTEX-A flights. Observed OH increases with altitude from 0.15 

15 pptv at the surface to 0.43 pptv at 1 1 km, while HO 2 decreases from 19 pptv at the surface to 5 

16 pptv at 1 1 km reflecting the positive relationship of OH/HO 2 ratio with NO [ Jaegle et al., 2000]. 

17 H 2 0 2 concentrations gradually decrease from 2.6 pptv at the surface to 0.15 pptv at 1 1km, 

18 reflecting the decrease in water vapor [Heikes et al., 1992], CTEO also decreases with altitude, 

19 reflecting the boundary layer source from oxidation of isoprene [ Millet et al., 2006], HNO4, a 

20 thermally unstable molecules that is an important HO x reservoir in the upper troposphere [Jaegle 

21 et al., 2000], has an 8-9 km peak of 62 pptv. 

22 OH and HO 2 in the original simulation match observations in the lowest 1 km, but are too 

23 high above by 60% and 30% respectively. Measurement nominal accuracy for both is 15%, well 

24 below this bias [Ren et al., 2006]. Increasing the lightning NO x source in the improved simulation 

25 decreases HO 2 while increasing OH. This fully corrects the HO 2 overestimate but worsens the OH 

26 overestimate above 5 km. Similar HO x simulation biases are found in a box photochemical model 
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1 constrained with local ICARTT observations [ Olson et al., 2006b; Ren et al., 2006], and which 

2 closely reproduces GEOS-Chem results. The OH bias must reflect either instrument error or a 

3 fundamental flaw in current understanding of HO x photochemistry; see Olson et al. [2006b] and 

4 Ren et al. [2006] for further discussion. In what follows we take the model OH at face value but 

5 will also discuss the effect of possible model bias. 

6 The model simulates in general well the concentrations of the HO x reservoirs H2O2, 

7 CH 2 0, and HNO4 (Figure 3). H2O2 is unbiased in the free troposphere, but is overestimated below 

8 3 km by ~30%; a possible cause is the model conversion of H0 2 to H 2 0 2 in aerosols [ Martin et 

9 al., 2003], which recent measurements find to be much lower than previously estimated 

10 [ Thornton and Abbatt, 2005]. The CH 2 0 simulation, discussed in detail by Millet et al., [2006], 

1 1 matches observations closely and is largely insensitive to changes in OH (which affect sources 

12 and sinks in the same direction). HNO4 is well matched in the original simulation but is 

13 overestimated by 30% when the lightning source is increased. 

14 5. Boundary Layer NO* over the United States 

15 Figure 4 compares simulated and observed mean vertical distributions of CO, NO x , PAN, 

16 HNO3, and ozone concentrations for the ensemble of DC-8 and WP- 3 D flights. Observed CO, 

17 NO x , PAN, and HNO3 are elevated below 3 km, reflecting anthropogenic sources. The model 

18 with original emissions is too high for all four species. Reducing the U.S. NO x emissions from 

19 stationary sources by 50% from 1999 to 2004 to account for the NO x SIP Call [ Frost et al, 2006] 

20 removes much of the model boundary layer NO x bias over the Midwest United States. This 

21 corresponds to a 20% overall in the anthropogenic NO x source (Table 1). Parrish, [2006], using 

22 urban ambient measurements of the CO/NO x ratio along with fuel sales data, previously found 

23 that the on-road vehicular source in the NE199 inventory is accurate for NO x but 50% too high for 

24 CO. Applying this correction to the on-road vehicular CO source, which represents 60% of total 

25 CO emissions in the United States according to NEI99, decreases the mean CO bias from 35 
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1 ppbv to 20 ppbv in the WP-3D data and from 25 to 19 ppbv for the DC -8 data (Figure 4). Further 

2 source reduction would be needed to match the boundary layer observations for CO. 

3 

4 6. Upper tropospheric NO, over the United States 

5 NO x concentrations observed in ICARTT show a decrease from the boundary layer to the free 

6 troposphere, but then a sharp rise with altitude above 6 km (Figure 4). Mean concentrations reach 

7 0.55 ± 0.36 ppbv at 8-12 km altitude, higher than in the boundary layer. The NO/NO x molar ratio 

8 averages 75% both in the observations and the model for the 8-12 km region (all data are for 

9 daytime). Observed PAN shows a broad maximum at 6-10 km. In contrast, HNO 3 is depleted in 

10 the free troposphere due to scavenging during uplift. Mean ozone increases with altitude from 50 

1 1 ppbv near the surface to 75 ppbv at 8 km. 

12 The original model greatly underestimates the upper tropospheric NO x enhancement; the 

1 3 discrepancy increases with altitude from a factor of 3 at 8 km to a factor of 5 at 1 1 km. Simulated 

14 PAN is too low by ~30% while HNO3 is well simulated. Ozone is too low by 10 ppbv throughout 

15 the free troposphere. 

16 Measurements of upper tropospheric NO x from previous aircraft campaigns over the United 

17 States indicate much lower concentrations than observed in ICARTT. Jaegle et al., [1998] report 

18 mean NO concentrations from the SUCCESS campaign out of Kansas (April - May 1996) of 

19 0.030 ± 0.022 ppbv for 8-10 km and 0.061 ± 0.045 ppbv for 10-12 km. Ridley et al., [1994 ] 

20 report mean NO concentrations of 0.2± 0.1 ppbv over New Mexico during during ELCHEM 

2 1 (July - August 1989), even though convection was frequently targeted. SONEX observations in 

22 October-November 1997 over Maine and Atlantic Canada indicate mean NO concentrations 

23 between 6-12 km of 0.1 ppbv (all data) and 0.23 ppbv (convective outflow) [ Crawford et al., 

24 2000]. SUCCESS and SONEX were in spring and fall, whereas lightning over the United States 
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1 peaks in June- August; ELCHEM took place in the southwestern United States whereas maximum 

2 lightning is over the Gulf Coast [ Orville and Huffmes, 2001]. 

3 Figure 5 (left) shows the observed spatial distribution of mean upper tropospheric NO x 

4 concentrations. Values exceeding 1 ppbv extend over much of the Southeast and Midwest. Deep 

5 convective injection of boundary layer pollution cannot explain these high values since the NO x 

6 mixing ratio above 8 km is greater than that in the boundary layer (Figure 4). Aircraft emissions 

7 cannot provide an explanation either because the geographical distribution does not match the 

8 aircraft corridor along the Eastern United States and North Atlantic [ Gauss et a!., 2006], and in 

9 any case these emissions are fairly well constrained from atmospheric measurements [Meijer et 

10 al., 2000] and represent only a small source of upper tropospheric NO x (Table 1). Transport of 

1 1 NO x from the stratosphere is also well constrained by atmospheric observations and is negligibly 

12 small, about 0.2 Tg N yr' 1 globally [Murphy et al., 1993; Jacob et al, 1996] 

13 Lightning provides the best explanation for the elevated NO x in the upper troposphere 

14 during ICARTT. Figure 6 (left) shows National Lighting Detection Network (NLDN) mean 

15 lightning flash rates for July 1- August 15. The NLDN data (>100 sites in the continental United 

16 States) were collected by Vaisala (www.vaisala.com) and supplied to us by the Global Hydrology 

17 Resource Center at NASA Marshall Space Flight Center. The NLDN network measures only 

1 8 cloud-to-ground lightning flashes, and intracloud flashes are estimated to be about 3 times that 

19 amount [Boccippio et al ., 2001]. Cloud-to-ground flash detection efficiencies are >90% over the 

20 continental United States and degrade rapidly off-shore and beyond U.S. borders. The NLDN 

2 1 lightning activity is heaviest in the Gulf of Mexico region but also has maxima in the Midwest. 

22 An additional tongue of lightning activity extends along the southwest U.S. over Arizona and 

23 New Mexico. We see substantial coincidence in Figure 5 and 6 between the geographical 

24 distribution of lightning and that of upper tropospheric NO x . 

25 2004 was not an anomalous year for lightning over the United States, as shown in Figure 

26 7 with total NLDN lightning flash counts for July-August 2000-2005. There is a notable jump in 
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1 lightning flashes from 2001 to 2002, which reflects an upgrade to detection completed during 

2 2002. The comparable years are thus 2003-2005. GEOS-Chem lightning counts computed using 

3 consistent GEOS-4 meteorology for 2000-2005 also show relatively little interannual variability 

4 in total lightning over the contiguous United States, with 2004 being typical. We find in the 

5 model that lightning flash rates over the eastern United States in 2004 are everywhere within 20% 

6 of the 1995-2005 mean. 

7 We can make a rough estimate of lightning emissions over the United States during 

8 ICARTT by using NLDN flash rates, which only measure cloud-to-ground, multiplying by 4 to 

9 account for intracloud flashes [Boccipio et al.„ 2001], and assuming a 500 mol NO x /flash 

10 production rate derived from the mean peak NLDN current [Ott et al., 2006]. We deduce an 

1 1 emission of 0.45 Tg N from the NLDN dataset. This is a factor of 7 above the standard GEOS- 

12 Chem simulation (Table 1). The model captures the maximum along the Gulf Coast but is too low 

13 offshore and over the Midwest (Figure 6). We tried to improve this model distribution with 

14 alternate lightning parameterizations based on cloud mass flux or convective precipitation [Allen 

15 and Pickering, 2002], but the cloud mass flux parameterization did not capture the Gulf 

1 6 maximum while the convective precipitation scheme did not capture the lightning distribution 

1 7 over land. 

18 The global lightning source of NO x in GEOS-Chem is 4.7 Tg N yr' 1 from 2.7 x 10 9 

19 flashes, corresponding to 125 mol/flash, a factor of 4 below the Ott et al., [2006] estimate, 

20 justifying a factor of 4 increase to the GEOS-Chem lightning source over the United States. 

2 1 Figure 5 (right) shows the resulting mean 8-12 km NO x concentrations, successful over the South 

22 where lightning flashes are correctly simulated, but still showing discrepancies in the Midwest 

23 due to model error in the geographical distribution of lightning. 

24 On a global scale though, a NO x yield of 500 mol/flash would appear to lead to an 

25 excessive lightning source. The OTD-LIS vl.O gridded satellite lightning climatology produced 

26 by the NASA L1S/OTD Science Team (Principal Investigator, Dr. Hugh J. Christian, NASA / 
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Marshall Space Flight Center), available from the Global Hydrology Resource Center 
(http://ghrc.msfc.nasa.gov) , yields 1.5 x 10 9 flashes yr' 1 . Combining a N0 X yield of 500 mol/flash 
with this global estimate would imply a lightning source of 10.5 Tg N/year, which seems too high 
based on constraints from satellite observations [ Boersma et al., 2005; Martin et al., 2006b] and 
tropical ozonesondes [ Martin et al., 2002]. While physical mechanisms responsible are not well 
understood, recent observational evidence suggests that NOx yields per flash are lower in tropical 
than in mid-latitude storms [Huntrieser et al, 2006]. 

The lifetime of NO x in the upper troposphere is a major uncertainty in scaling the 
lightning source to match the ICARTT observations. As shown in section 4, the model OH 
concentration in the upper troposphere is a factor of 2 higher than observed. Figure 8 shows 
simulated and observed frequency distributions of NO x concentrations at 8-12 km on the 2°x2.5° 
model grid. If the model lifetime were too short due to excessive OH, then one might expect the 
variability in the model to be larger than observed, but Figure 8 shows that this is not the case, 
adding some support to the simulated NO x lifetime. 

The large lightning source inferred from the ICARTT observations has important 
implications for tropospheric ozone. Li et al., [2005] found that a semi-permanent upper level 
cyclone above the southern United States in summer allows ozone build-up in the upper 
troposphere by trapping convectively lifted precursors and lightning NO x . Cooper et al., [2006] 
confirmed the resulting ozone maximum by analysis of ozonesonde data during ICARTT, and 
found an associated 11-14 ppbv contribution to ozone from lightning averaged over the eastern 
United States. Figure 4 shows that the improved simulation with increased lightning largely 
removes the upper tropospheric ozone bias in the original simulation (reduction in surface 
anthropogenic NO x emissions in the improved simulation decreases upper tropospheric ozone by 
only 1-2 ppbv). The residual bias appears due to insufficient lighting generation in the Midwest 
(Figure 6). We find in the model that lightning enhances upper tropospheric ozone by ~10 ppbv, 
consistent with the analysis of Cooper et al., [2006]. 


15/29 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


7. Chemical evolution and export of U.S. NO x emissions 

In this section we use the ICARTT data to estimate the export of anthropogenic NO y from the 
U.S. boundary layer to the free troposphere and the speciation of this NO y . Figure 9 shows the 
simulated and observed NO y speciation below 2km. Here and from now on model results are 
from the improved simulation with reduced fuel NO x emissions and increased lightning (Table 1). 
The mean observed NO x /NO y ratio at 0-2 km is 30% in the Northeast/Midwest, reflecting the 
density of sources, 21% in the South, and 15% offshore. The dominant component of NO y in all 
three regions is HNO 3 , averaging 53% in the Northeast/Midwest and 75% offshore. PAN is 
favored in the South, despite high temperatures, likely due to large isoprene emission [ Horowitz 
et al., 1998]. The mean PAN/NOx ratio is 1.02 in the Northeast/Midwest and 1.54 in the South. 
The model is remarkably successful at reproducing these fractions and patterns. 

We can estimate the export efficiency / of NO y from the North American boundary layer, 
following the approach of Parrish et al. [2004], by viewing CO as an inert tracer and comparing 
the enhancement ratio ANO y /ACO in North American pollution outflow to the anthropogenic 
molar emission ratio R of CO to NO x : 


f = R-a 


AM?y 

A CO 


( 1 ) 


where a is a correction factor to account for the CO source from boundary layer 
oxidation of nonmethane hydrocarbons, particularly biogenic isoprene [ Chin et al., 1994]. 
Derivation of the export efficiency following (1) requires estimates of R , a, and the 
background NO y and CO concentrations to which the enhancements A are referenced. By 
adopting constant values for these variables, as discussed by Li et al. [2004], we obtain a 
simple observationally-based diagnostic of export for testing more elaborate models such 
as GEOS-Chem. We use here and background concentrations of 95 ppbv CO [Li et al., 2004] 
and 100 pptv NO y [Parrish et al., 2004], a = 1.2 [Chin et al., 1994], and R = 5.9 mol mol ' 1 from 
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1 our improved GEOS-Chem simulation for the United States east of 100°W. R varies little by 

2 region between the Midwest (5.6), South(5.9), and Northeast(6. 1). 

3 We apply equation (1) to every anthropogenic pollution plume observed by the DC-8 and 

4 WP-3D between 2.5 and 6.5km as defined by a CO enhancement ACO > 30 ppbv. Plumes above 

5 6 km are ignored due to lightning NO x interference, and we also exclude biomass burning plumes 

6 diagnosed from nitrile data (section 2). Anthropogenic pollution plumes defined in this manner 

7 represent 1 1% of the combined INTEX-A and ITCT2k4 data at 2. 5-6. 5 km. From these data we 

8 find a mean NO y export efficiency f =16±10% to the free troposphere with mean composition of 

9 13% NO x , 40% PAN, and 47% UNO,. Sampling the model along the ICARTT flight tracks 

10 shows a comparable value in both magnitude and variability: /=14 ± 8% with mean composition 

11 of 9% NO x , 42% PAN, and 49% HNO3. Previous studies using aircraft data for North American 

12 outflow in NARE’97 [Li et al., 2004; Parrish et al., 2004] and Asian outflow in TRACE-P 

13 [Koike et al., 2003; Miyazaki et al., 2003] similarly found / values in the range 10-20%. 

14 Figure 10 shows the vertical distribution of / and the speciation of this exported NO y . Highest 

15 observed mean and variability of/(18 ±11%) is at 2. 5-3. 5 km, where HNO3 dominates the NO y 

16 export fraction (54%). The model in that altitude range shows a similar value of/(15 ± 9 %) and 

17 HNO3 fraction (55%). Parrish et al. [2004] proposed that this shallow venting is due to fair 

1 8 weather cumulus breaking through the afternoon boundary layer. In the model, outflow from such 

19 shallow wet convection would experience only limited scavenging of soluble species [Liu et al., 

20 2001]. In a 2-km deep updraft column 63% of HN0 3 would be scavenged, allowing for some 

2 1 HNO3 export to the free troposphere, consistent with observations. 

22 The export efficiency/ decreases with altitude in the model and observations 

23 reflecting HNO 3 scavenging during uplift, except between 5. 5-6. 5 km in the model, 

24 where lightning may begin to affect the calculation. We see from Figure 10 that PAN is the 

25 principal component of exported NO y above 4 km, both in the observations and in the model. 
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1 This dominance of PAN in free tropospheric continental outflow of NO y has been previously 

2 observed in aircraft campaigns downwind of North America [. Parrish et al., 2004] and Asia 

3 [Miyazaki et al., 2003], 

4 8 Conclusions 

5 The ICARTT study in summer 2004 provided extensive observations of reactive nitrogen 

6 (NOy) species over the eastern United States and western North Atlantic, from the surface to 12 

7 km altitude. We have interpreted these observations with a global 3-D model of tropospheric 

8 chemistry (GEOS-Chem) to place constraints on the sources, chemical evolution, and export of 

9 NOy from North America. 

10 ICARTT observations in the continental boundary layer provide top-down verification of the 

1 1 recent decrease in stationary NO x emissions in the eastern United States mandated by the NO x 

12 SIP Call. Model simulation of NO x in ICARTT indicates that the latest comprehensive national 

13 emission inventory done for 1999 (NEI 99) is too high over the Midwest by almost a factor of 2. 

14 This is consistent with the 50% reduction in stationary sources from 1999 to 2004 inferred from 

15 power plant smokestack monitoring [ Frost et al., 2006] and amounts to a 22% decrease in U.S. 

16 anthropogenic NO x emissions, to 0.62 Tg N for the July 1-August 15 ICARTT period. GEOS- 

17 Chem emissions of anthropogenic NO x during that same period were 0.98 Tg N for East Asia and 

18 0.53 Tg N for Europe. 

19 Observed NO x concentrations in ICARTT show a sharp rise above 6 km with mean 

20 concentrations reaching 0.55 ± 0.36 ppbv at 8-12 km, higher than observed in the U.S. boundary 

2 1 layer and much higher than observed in the upper troposphere on previous U.S. aircraft 

22 campaigns (SONEX, SUCCESS, ELCHEM). A close correspondence is observed between the 

23 spatial distribution of upper tropospheric NO x during ICARTT and lightning flash counts 

24 observed by the National Lightning Detection Network (NLDN), identifying lightning as the 

25 dominant source for the observed NO x . Using NLDN flash rates for July 1- August 15, scaled up 

26 by a factor of 4 to account for intra-cloud flashes [Boccippio et al., 2001], and assuming a 500 
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1 mol NO x /flash production rate following Ott et al., [2006], we deduce a lightning NO x emission 

2 of 0.45 Tg N over the United States and adjacent coastal areas. This is a factor of 7 higher than in 

3 the standard GEOS-Chem simulation, which uses a yield of 125 mol NO x /flash for a global 

4 lightning NO x source of 4.7 Tg N yr" 1 . We can reproduce the upper tropospheric NO x observed 

5 over the southern United States in ICARTT by increasing the lightning NO x yield in the model by 

6 a factor of four to the Ott et al. [2006] value. 

7 This also provides a successful simulation of the observed frequency distribution of upper 

8 tropospheric NO x , and corrects a 5-10 ppbv low bias in the simulation of ozone in the free 

9 troposphere in the model. The factor of four increase is probably not extrapolatable globally as 

10 the resulting global lightning source in the model (19 Tg N yr' 1 ) would be too high relative to 

1 1 observational constraints [ Boersma et al., 2005; Martin et al., 2006b], Recent observations 

12 suggest that mid-latitude storms have higher NOx yields than tropical storms [ Huntrieser et al., 

13 2006]. 

14 Uncertainty in OH concentrations in the upper troposphere is a significant limitation for 

15 interpreting quantitatively the observed NO x concentrations in terms of an implied lightning NO x 

16 source. Simulated upper tropospheric OH concentrations in GEOS-Chem are about a factor of 2 

1 7 higher than observed in ICARTT, and the same bias is found in box photochemical model 

18 calculations constrained with the aircraft observations [ Olson et al., 2006b; Ren et al., 2006]. As 

19 discussed by Olson et al. [2006b] and also in Spivakovsk y et al. [2000], this discrepancy is 

20 beyond what one might expect from standard error propagation in a photochemical model, and 

2 1 points either to instrument error or to some fundamental flaw in understanding of upper 

22 tropospheric HO x chemistry that would in turn affect the simulated NO x lifetime. A model 

23 decrease of OH by a factor of 2 would correspondingly decrease the required increase in the 

24 lightning NO x source. 

25 We examined the speciation of NO y over the United States, in the observations and in the 

26 model, to gain insight into the chemical evolution and export of NO y . The NO x fraction is largest 
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1 in the boundary layer over the Northeast, averaging 30%, and reflecting the density of sources. 

2 PAN makes a somewhat larger contribution to NO y in the Southeast than elsewhere, possibly 

3 reflecting isoprene emissions. The model reproduces well the observed partitioning of NO y for all 

4 regional and altitudes, implying a good understanding of NO y chemistry. 

5 We estimated the export efficiency / and related speciation of NO y out of the North American 

6 boundary layer with a Lagrangian analysis of NO y -CO correlations in the free troposphere, 

7 following the approach of Parrish et al. [2004], For the ICARTT data at 2. 5-6. 5 km altitude we 

8 find /= 16±10% in the observations and 14±8% in the model, consistent with previous studies of 

9 North American and Asian outflow [Li et al., 2004; Parrish et al., 2004; Koike et al., 2003; 

10 Miyazaki et al., 2003]. The highest export efficiency is in the lower free troposphere but is then 

1 1 mostly HNO3 venting from shallow convection. We find that PAN is the dominant component of 

12 exported NO y (>50%) in pollution plumes above 3.5 km, consistent with previous studies of 

13 Asian outflow. The successful simulation of export of North American NO y offers confidence in 

14 current model estimates of North American influence on the global NO y and ozone budgets [Li et 

15 al., 2004] with implications for intercontinental pollution [Li et al., 2002].. 

16 
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1 Table Caption 

2 

3 Table 1. Contiguous United States NOx Emissions for July 1 - August 15, 2004 

4 

5 Figure Captions 

6 

7 Figure 1. Flight tracks for ITCT 2k4 (right) and INTEX-NA (left) aircraft campaigns 

8 (July 1- August 15) shaded by altitude. Shaded areas define regions (left) used in 

9 comparisons between model and observations: South (mid-grey), Northeast (light grey) 

10 and Midwest (dark grey). 

11 

12 Figure 2. Mean July-August anthropogenic NO x emissions in the United States. 1999 

13 values from the EPA National Emissions Inventory (NEI99 vl) are compared to 2004 

14 values including a 50% reduction in point source emissions. Emissions from Canada and 

15 Mexico are from GEIA scaled to 1998 following Bev et al., [2001]. Colorscale saturates 

16 at 5x10" molecules/cm" s’ 

17 

18 Figure 3. Mean vertical profiles of OH, HO2, H2O2, HCHO and HNO4 concentrations. 

19 Observations from the DC-8 aircraft (thick black) are compared to model results using 

20 the original (dashed grey) and improved (thick grey) simulations and compared to 

21 observations (black). Horizontal bars are standard deviations on the observations. 

22 Here and in subsequent figures, the ICARTT observations have been filtered to remove 

23 urban plumes, biomass burning plumes, and stratospheric air as described in text. Model 

24 results are sampled along the flight tracks at the time of flights. Modifications to 

25 emissions from the original to the improved simulation include a four-fold increase in the 

26 lightning source, upward extension of the lightning source to the July mean tropopause 

27 height, and a 50% reduction of the NEI99 vl CO transport and NO x point source 

28 emissions. 

29 

30 Figure 4. Same as Figure 3 for CO, NO x (median), PAN, HN03,and ozone 

3 1 concentrations as sampled by the DC-8 (top) and the WP-3D (bottom) aircraft. 

32 

33 Figure 5. Mean upper tropospheric NO x concentrations (8-12 km) during ICARTT (July - 

34 August 2004). Observations (left) are compared to improved model values (right) over 

35 the 2ox2.5o model grid. The improved model has a factor of four increase in the United 

36 States lightning emissions relative to the standard GEOS-Chem version. 

37 

38 Figure 6. Mean lightning flash rates for July 1 - August 15, 2004. Observations of cloud- 

39 to-ground lightning from the National Lightning Detection Network, multiplied 

40 uniformly by a factor of 4 to account for intracloud lightning, are compared to standard 

41 GEOS-Chem results. 

42 

43 Figure 7. Cloud-to-ground lightning flash counts in July-August 2000-2005 over the 

44 United States. National Lightning Detection Network flash counts (grey bars) are 

45 compared to modeled flash counts in GEOS-Chem derived using GEOS-4 meteorology 

46 for the domain 130-70°W, 25-50°N (grey and white bars). Model flash counts are divided 
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1 uniformly by a factor of 4 to account for intracloud lightning [Boccipio et al., 2001]. The 

2 jump in NLDN data between 2001 to 2002 reflects an upgrade in detection; thus, 

3 comparable years are 2003-2005. 

4 

5 Figure 8. Frequency distribution of concentrations at 8-12 km altitude during ICARTT. 

6 Observations (line) are compared to results from the improved model with increased 

7 lightning source (grey bars). The frequency distribution is shown as he number of 

8 occurances on the 2°x2.5° grid. 

9 

10 Figure 9. Mean relative molar contribution of NO x , PAN and HNO3, to total NO y defined 

11 as NOy=NO x +PAN+HN03 in the ICARTT data (July 1 - August 15, 2004). Urban, 

12 biomass burning, and stratospheric plumes have been excluded as described in the text. 

13 Observations from the DC-8 and WP-3D aircraft are compared to model results 

14 sampled along the flight tracks. 

15 

16 Figure 10. Speciated export efficiency of NOy (defined as NO y = NO x + PAN + HNO3) 

17 from the North American boundary layer during ICARTT, derived from equation (1) as a 

18 function of altitude. Observations (left) are compared to model results (right). 

19 



